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Abstract 

A numerical model has been developed which can predict 
both the unsteady flows within a wave rotor and the 
steady averaged flows in the ports. The model is based on 
the assumptions of one-dimensional, unsteady, perfect gas 
flow. Besides the dominant wave behavior, it is also 
capable of predicting the effects of finite tube opening 
time, leakage from the tube ends, and viscosity. The 
relative simplicity of the model makes it useful for design 
and optimization, as well as analysis, of wave rotor cycles 
for any application. This paper will discuss some details 
of the model and will then present comparisons between 
the model and two laboratory wave rotor experiments. 

I .0 Background 

The Wave Rotor is presently being investigated for use as 
a core gas generator for multi-spool gas turbine engines in 
order to achieve very high peak cycle temperatures and 
pressures. The device, shown schematically in Figure 1.1, 
uses unsteady waves to transfer energy directly to and 
from the working fluid through which the waves travel. 
It consists of a series of constant area passages (tubes) that 
rotate about an axis. Through rotation, the ends of the 
tubes are periodically exposed to various circumferentially 
arranged ports in which the flow is steady, but which 
initiate unsteady waves within the tubes. Because of the 
unsteady nature of the device, each tube of the wave rotor 
is periodically exposed to both hot and cold flow over 
roughly equal time periods. The mean temperature of the 
rotor material may therefore be expected to remain 
considerably below the peak cycle temperature. This 
characteristic, a comparatively low rotational speed for a 
given mass flow (as determined through analysis), and the 
absence of high torque loads, give the wave rotor potential 
advantages over the conventional gas generator. 

Analyzing the wave rotor is difficult for several reasons. 
First, although the flows in the ports are steady, there is 
no steady state analysis which allows their properties to be 
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determined. The wave rotor must therefore be analyzed 
by following the unsteady fluid dynamics that occur within 
one or several of the passages (tubes) as they rotate past 
the ports and walls. Second, the governing differential 
equations of motion in the tubes are hyperbolic and often 
have shock waves present. Closed form solutions of these 
equations do not generally exist, so numerical methods 
must be used. An effort is presently underway by the first 
author to develop such a numerical model or code. The 
initial phases of this effort were reported previously 0 *. In 
that report the initial simplifying assumptions for the 
model were presented. For reference these were uniform 
cross sectional (i.e. one-dimensional), inviscid, perfect gas 
flow, with tubes which were instantaneously exposed to 
the ports in the manner of a ruptured shock tube 
diaphragm. Furthermore, it was assumed that since each 
tube in the wave rotor undergoes the same cycle, 
following one tube was equivalent to following them all. 
This amounts to an assumption that the tubes on the rotor 
do not interact. These simplifications reduced the 
modelling of the wave rotor to integrating the Euler 
Equations. A comparison was made between the model 
and an actual experimental wave rotor. It was found that 
although the model was qualitatively correct, there were 
significant discrepancies from the experiment. 

The present paper is intended to highlight some of the 
modifications that have been made to the original model 




in order to predict the wave rotor behavior more 
accurately. In particular, the model, though still one 
dimensional, is now capable of including the effects of 
viscosity, leakage from the ends of the wave rotor tubes, 
and non-ins tantaneous opening of the tubes to the ports. 
These effects have been studied to some degree by other 
investigators. In particular, Eidelman' 2 ’ studied finite 
opening time effects using a two dimensional inviscid 
code, and Taussig'-'' investigated all three effects with a 
one dimensional model. However, Eidel man’s code is 
computationally expensive and thus not practical for 
analyzing entire wave rotor cycles. Furthermore, it does 
not account for leakage or viscosity. Taussig’s code is not 
described in the literature with enough detail to 
implement. Furthermore, he assumes that the leakage 
effects are uncoupled from the wave phenomena, which 
from the present study, does not appear to be the case. 
The paper will proceed by discussing the three 
modifications outlined above in some detail, then 
comparing the modified model to two simple wave rotor 
experiments. 


2.0 Finite Opening Time 

It is often assumed in wave rotor analysis that the ends of 
the tubes are instantaneously opened or closed as they 
rotate into or out of the various ports. This may be 
interpreted as assuming that the time required for the tube 
to pass from a fully open to a fully closed position (or 
vice-versa) is much less than the time required for a 
characteristic wave to travel down the length. In practice 
this is almost never the case, and the effects on the wave 
rotor performance due to this so called finite opening time 
can be substantial. Of course, the finite opening time 
effects are highly three dimensional; however, it is 
possible to obtain a reasonable estimate of them by 
modifying the boundary conditions for the one 
dimensional model. The discussion to follow will 
consider the left end of the wave rotor tube; however, it 
is understood that analogous conditions exist for the right 
end. 

Consider Figure 2.1, where the left end of the tube is 
shown as it begins to enter the port region. Depending 
upon the port conditions, the flow will be either in or out. 
Looking first at the inflow scenario (a) it is imagined that 
the flow in the port itself is isentropic relative to some 
specified stagnation state, and that it instantaneously 
adjusts to conditions downstream (i.e. steady nozzle flow). 
Adjacent to this region is an imaginary "mixing region” 
where the fluid changes from a highly non-uniform 
distribution across the partially open tube on the left to 
some mixed out average on the right. The mixing is also 
imagined to occur instantaneously, and the extent of the 
mixing region is assumed negligible. Adjacent to this 
region is the first "cell” of the computing domain for 
which, at the present time, the conditions are completely 
known. Assuming that the flow is fiictionless in the 
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Figure 2.1 Finite Opening Time Model 


mixing zone, and that the pressure is uniform across the 
left face, the following mass, momentum, and energy 
equations may be written for a perfect gas: 

P« u o = P e « e ^ * 4> m (2.1) 

, , A, 

Po + Po u o = P e + PX-^ = <f> m „ ( 2 . 2 ) 

^I'uPo+^Po^ = Pe u o H -^ - <tV ( 2 - 3 ) 

where H is the total enthalpy. The subscripts 0 and e in 
these equations designate the end and the beginning planes 
of the mixing zone respectively, as shown in Figure 2.1. 
If the pressure at the exit of the isentropic region, p e is 
specified, then given the known stagnation conditions of 

A 

the duct and the area ratio of the exit to the tube, I, the 

A ' 

conditions at the exit are completely known and the 
quantities (|) J0 , <J> mo , and (|> c in equations 2. 1-2.3 can be 
calculated. Equations 2. 1-2.3 may also be combined to 
yield the quadratic 
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This may be solved for (choosing the subsonic value). 
With u 0 known, equations 2.1 and 2.2 may be used to 
obtain p () and p 0 so (hat the state at the right end of the 
mixing zone is completely known. 

The conditions p„ p,, and u, in the first computational cell 
of the tube are also known. The difference between p 0 
and p, will give rise to either a compression or expansion 
wave travelling to the right. The pressure ratio across the 
wave and the conditions in the first cell are sufficient 
information to analytically determine the velocity behind 
it, u 0 \ using either the normal shock relations or isentropic 
centered expansion equations' 4 *. Since nothing has been 
said about the choice of p 0 it is not expected that u 0 * will 
be the same as calculated using equations 2. 1 -2.3. A 
function may be defined however as 

y(Pc) =u o- u o‘- (2.5) 

The proper choice of p e is that for which y=0. This 
cannot be found analytically but may be obtained using a 
convenient numerical root finding technique (e.g. the false 
point method (S) ). The values of p,„ p () , and ii*, obtained 
through the solution of equation 2.5 are then assigned to 
the left image cell of the computing space. 

It is possible that the pressure p c drops to the point where 
the flow at the exit of the isentropic region is sonic but 
where equation 2.5 has not yet reached zero. In this case, 
it can no longer be assumed that is uniform across the 
tube and equation 2.2 can no longer be used, Equations 
2.1 and 2.3 remain valid however, and furthermore the 
values of <!>,„ and <|\. become frozen at their choked limits. 
As such, equation 2.3 may be rewritten as 
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Thus, for any value of p 0 chosen, u 0 is known. Again, the 
conditions in the first cell are known so that another 
function 


y(Po) =su ir u <>* (2.7) 

may be defined and the value of p 0 for which y=0 may be 
found through the same method as equation 2.5. 

For outflow, consider Figure 2.1(b). Here, the duct is 
assumed to be a constant pressure region with the 
pressure, p c known. Adjacent to this is a fictitious zone in 
which the flow is isentropic and which, like the mixing 
zone above, adjusts instantaneously to flow conditions. 
Again, conditions in the first cell of the computing space 
are known. Assuming that the flow is unchoked at the 
exit of the isentropic zone, and subsonic in the first cell, 


steady isentropic flow laws require that between planes 0 
and e, the Mach numbers are related by 
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Thus, if p 0 is specified then the exit Mach number is 
known. Since the state in the first cell is completely 
known, specification of p 0 also completely determines the 
state at the right plane of the isentropic zone (either 
through normal shock laws or laws for centered fans as in 
the inflow problem above). With p () , u 0 , and p 0 known, 
along with p tf and the exit Mach number, M t „ u e and p e 
may be found using the relations for steady isentropic 
flow. Again, nothing has been said about the choice of p 0 ; 
however, in keeping with the assumptions laid down, the 
correct value is that for which mass continuity is 
maintained across the isentropic zone, namely 


y(P„) = Pc U e-^-Po ll O * 0 (2.9) 

As with the inflow problem, equation 2.9 must be solved 
numerically using some iterative method. The resulting 
values of the primitive variables for plane 0 are then 
assigned to the image cell of the computing space. If, in 
the process of solving equation 2.9, the exit Mach number 
drops below a value of -1.0 then the isentropic region is 
assumed choked and tire values of p e and u e corresponding 
to this condition are used. 

Additional possibilities must be considered for the outflow 
condition if the flow is negatively supersonic in the first 
cell. In this case p 0 is assigned a value, p {) * which would 
create a stationary shock in plane O' 4 *. The state at the 
exit is then found using equation 2.8 and the steady 
isentropic flow laws. The value of y in equation 2.9 is 
then calculated, if this value is greater than zero then the 
iteration process may continue since it is now assured that 
there is some p e < p 0 < p 0 * t for which ;y-0. If y is less than 
zero the only remaining possibility is that the flow in the 
first cell is completely unaffected by the boundaries. In 
this case, the image cell will have no effect on the internal 
computations and it may simply be assigned the values of 
the first cell. 


In the above discussion it is assumed known a priori 
whether the flow resulting from the imposition of 
boundary conditions will be inwaid or outward. In 
actuality, the code must anticipate the resulting flow and 
then apply the appropriate model. This may be done in a 
rather straightforward manner as follows. If the flow is 
assumed inward, then the minimum value of equation 2.5 
is obtained by setting p e equal to the prescribed duct 
stagnation pressure, in this extreme case, p e is also the 
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same as p 0 and the value of y is not only minimum, but 
negative. The exit pressure, p e must be lowered from this 
value in order to bring equation 2.5 to zero. If, the value 
of y is found to be positive, then the assumption of inflow 
was incorrect, and the outflow model must be used. It is 
noted here that in the code, the same prescribed boundary 
pressure is used whether the flow is in or out. Thus, for 
inflow, the pressure is interpreted as a stagnation value 
and an additional stagnation temperature is needed. For 
outflow, it is interpreted as a static pressure. An analogy 
may be made here if the wave rotor duct is viewed as a 
large reservoir or "tank". Flow into this tank simply loses 
its kinetic energy at the tank stagnation pressure (i.e. jet 
flow). Flow from the tank leaves isentropically at the 
stagnation pressure. 




pressure in the tube and the time has been normalized by 
the quantity L/a, where L is the tube length and a is the 
initial speed of sound in the tube. For this calculation a 
cell spacing of Ax/L=.02 was used, with an associated 
time step of Ata/L=.008. Figure 2(a) illustrates the 
scenario with instantaneous opening time. JHere, a shock 
wave is formed immediately and travels rightward with a 
uniform speed u sh =1.65. Figure 2(b) shows the same 
scenario with an opening time t O|)eu a/L=0.3. Although it is 
not possible to claim quantitatively that the results are 
correct, it is fair to say that the applied boundary 
conditions appear to generate the anticipated area-averaged 
behavior in the tube. In particular, there is a certain 
amount of time required for an actual shock to form, and 
the speed at which it travels is slower than that in the case 
where the opening was instantaneous. It is this effect on 
wave speed that is critical for wave rotor design and 
analysis and, as will be shown later, the approach 
described above produces excellent agreement with 
experimental data. 

3.0 Source Terms 

Aside from the effects of finite opening time, the fluid 
mechanics in the wave rotor tube are dominated by one 
dimensional, inviscid wave phenomena. As such, the 
Euler equations serve as a good starting point for 
predicting the time dependent behavior; however, there are 
several effects which, while secondary, are not negligible. 
In particular, it has been found that viscous effects, and 
the effects due to leakage of flow from the ends of the 
wave rotor tube to the casing can have a significant impact 
on the behavior of the flow. As with finite opening time, 
these effects are three dimensional; however, the one 
dimensional model can again be modified to account for 
them satisfactorily. The approach for both phenomena is 
to add source terms to the original homogeneous 
governing equations. The discussion to follow will first 
consider the addition of source terms in general to the 
Euler equations, and the requirements for high order 
accuracy. Following this, the particular form of the source 
terms used for modelling viscous effects (friction) and 
leakage effects will be presented. 

The particular form of the modified Euler equations to be 
considered is one in which the source term, like the flux 
vector, is a function of the original conserved vector. 
That is, 
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Figure 2.2 Finite Opening Time Effects 
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(3.1) 


As an example of the application of finite opening time 
boundary conditions, Figure 2.2 shows a computed time 
dependent pressure distribution in a hypothetical wave 
rotor tube as the left side is exposed to a high pressure 
duct. The pressure has been normalized by the initial 


For reference, the vectors w and F have the respective 
perfect gas forms: 
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equation 3.7 may be rewritten as 
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where h and e are the specific enthalpy and internal 
energy respectively, and the numbered subscripts refer to 
the elements of w. The form of S(w) has yet to be 
specified. 
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The straightforward approach to numerically integrating 
these equations is to expand all of the derivatives about a 
point, i.e. 


Multiplying equation 3.10 by At/2, replacing the spatial 
derivatives on the right with their numerical equivalents 
(note that these need only be first order accurate), and 
adding this to equation 3.6 yields 
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where the subscript n indicates the discrete temporal index 
nAt and i the spatial index iAx. Adding equations 3.4 and 
3.5, then subtracting the source term S, from both sides of 
the resulting equality yields 


w"* 1 - w" F“ -F” aw" aF" 

—I — i + --**' — s = — 1 + zl -s 
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The numerical scheme on the left hand side approximates 
the actual differential equation on the right up to a first 
order truncation error of order O(At). This can be 

0 2 w u 

improved by finding a suitable approximation for _Z!_. 

dt 2 

After Lax and Wendroff 6) , taking the time derivative of 
equation 3.1 at the point nAx, iAt yields 
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at 2 axat at v ' 

where the subscripts have temporarily been dropped. 
Using the Jacobians defined as 
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Here, the Jacobian [A] l+J/2 is the mean of two adjacent 
cells, i.e. %([A], +j +[A]j). The left hand side of this 
equality, when set equal to zero is the second order 
accurate numerical approximation to equation 3.1. This 
may be reairanged and compacted to the form 


w -'*' = w "-(f" - f» + s »At (3.12) 

— i — j \~4+i/2 -i-i aj Ax "~ l 
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(3.14) 


The numerical flux estimate for the cell face, fj +l/2 in 
equation 3.13, has been chosen in this particular fashion 
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because it is identical to the conservative form found when 
numerically approximating the homogeneous version of 
equation 3.1. It is found in practice however, that this 
estimation produces physically unrealistic oscillations in 
the solution. To overcome this problem, equation 3.13 is 
replaced with a high resolution, upwinding flux estimate 
guaranteed to produce monotonic results in the 
homogeneous Euler equations. This is based upon Roe’s (7) 
approximate Riemann solver and is described in detail in 
a previous report (1) . The numerical source estimate, ^ in 
equation 3.14 is seen to be based entirely on central 
differences. This is physically intuitive since the sources 
should not have any directional bias. There is no way to 
formally guarantee that the resulting scheme is monotonic 
since equation 3.1 is not conservative (see Lax [3]); 
however, the approach taken above will at least approach 
monotonicity as the strength of the source term is brought 
to zero. Furthermore, it is found in practice that this 
method retains very sharp resolution of discontinuities. 
This is illustrated in Figure 3.1 where the method has been 
applied to the scalar equation 


du 

~di 


+ u_ 


du 

dx 


= cu |u|. 


(3.15) 


The figure shows the initial profile, a half sine wave, 
corresponding to t=0. Also shown in the figure are the 
profiles at t=0.5 for three different values of the constant 
c in equation 3.15, corresponding to positive, negative and 
zero source term strengths. For this computation spatial 
and temporal spacings of Ax=0.0! and At-0.005 
respectively were used. There is no analytical solution by 
which to verify the results; however, it is seen that the 
shock remains sharp, and the pulse height and speed vary 
as expected. Furthermore, although not shown here, the 
stability of the scheme appears unaffected since the 
maximum allowable CFL number was found to be 
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Figure 3.1 Burger’s Equation With a Source 


3.1 Viscous Source Terms 


The major effect of viscosity in the wave rotor tube is that 
it creates boundary layers near the walls. These in turn 
alter the behavior of the core flow which remains inviscid. 
Naturally, a one dimensional model cannot analyze this 
phenomenon per se; however, if the terms in the one 
dimensional equations of motion are viewed as integral 
averages of the values over the tube cross section, then the 
viscous effects (except in the vicinity of shock waves) 
become confined to the shear stress at the wall, and the 
momentum equation takes the form 

4-pu+4-(p + p u2 ) = (3.i.i) 

at ax D. 


Here, D, is the hydraulic diameter defined as where 

W 

V 

A is the cross sectional area of the tube and W is the 
wetted perimeter. The task now becomes one of obtaining 
a reasonable estimate of Some researchers [4] have 
used correlations corresponding to fully developed 
turbulent pipe flow; however, considering typical tube 
dimensions and the amount of time that the gas in the tube 
is actually moving, this does not seem to be an appropriate 
approach. An alternate estimate of the wall shear stress 
may obtained by considering each discreet computational 
region of the tube to be a suddenly decelerated flat plate 
(or inversely a stationary plate with suddenly accelerated 
fluid). The solution to such a problem, for incompressible 
flow, is well known (8) and has the form 

u = U^erfitn) 

n = y - (3.1.2) 

2 v/vT 

Here, LL is the speed of the freestream flow, v is the 
kinematic viscosity, y is the distance from the wall, and 
erf is the error function. For the problem at hand, the 
freestream velocity is replaced by the local velocity in the 
tube. The wall shear stress has the form 
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(3.1.3) 


A characteristic time associated with the wave rotor is that 
required for a wave to travel down the tube, namely t=L/a* 
where L is the tube length and a' is the speed of sound at 
some characteristic stagnation state. Substituting this into 
equation 3.1.3, assuming constant viscosity, and non- 
dimensionalizing by p*, p*, and a\ the momentum 
equation may now be written as 
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Since the above analysis is only an estimate of the wall 
shear behavior, the term oc appearing in equation 3.1.4 has 
been added as a constant to be determined through 


comparison with experimental data. The ratio 


v,L 


D u 


is 


a measure of the penetration depth of the boundary layer 
into the tube and is generally small. When multiplied by 
the additional constant terms however, the contribution to 
the right hand side of equation 3.1.4 becomes significant 

(on the order of 0.1). 


3.2 Leakage Source Terms 

Leakage in the wave rotor occurs from the tube ends into 
the rotor casing and vice versa, as illustrated in Figure 
3.2.1. As such, the leakage source terms are applied only 
to the first and last cells in the computing domain and 
they affect only the continuity and energy equations. 
Consider the control volume shown in Figure 3.2.2. As 
with the boundary conditions, it is assumed that the 
leakage flow adjusts instantaneously to the conditions in 
the cell. The leakage gap is imagined as distributed 
uniformly along the cell in the manner of porosity. 
Assuming that the flow is out of the tube, the isentropic 
nozzle flow equation may be used to write 
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where is the leakage area per unit length of the tube, 

dx 

and p rav is the pressure in the center cavity of the wave 
rotor. Incorporating this equation into the continuity and 
energy equations yields 


dw, dF. 

+ i 

dt dx 
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where C, is the leakage coefficient C, = 


2y 


i/2 
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dx 


and A is the cross sectional area of the tube. The value 

f 
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( 
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of 


dx 


may be estimated by assuming a rectangular 


tube cross section in which case it takes on the value 


25 

, where 5 is the leakage gap, h is the tube height, and 

hAx 

Ax is the numerical cell size. If the flow becomes choked 
in the nozzle equation then the value of the flow function 
f in equation 3.2.1 becomes frozen at the choked limit. If 
the cavity pressure is higher than the cell pressure, then 
the pressures and densities in equation 3.2.1 are replaced 
with the cavity values and the cavity pressure in the flow 
function is replaced with the cell pressure. 



[ 

□mnrnrzz j 

p« n i 

X P \ 

X * 

T 

r 

r r 

leakage How 

r 

T 

P 1)11 XfAx 

X+Ax 

l >M x 1 






pll 

x+'ix 



"""x j 


Jl 

..k.J... 

X. 

X... 

pull 

X + AX 

□r r 


Ax 


Figure 3.2.2 Leakage Flow Model 
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It is not obvious from equations 3.2.2 and 3.3.3 that the 
source term is a function of the conserved vector w. More 
importantly, the form of the Jacobian [B] is difficult to 
obtain directly. It can however, be obtained through use 
of the chain rule, namely 

[ B ] = = a| an 

5w 5u <3w 

where u is the primitive variable vector (p, p, u) T . 

4.0 Results 

The results to be presented will compare predictions from 
the model described in the preceding sections with two 
laboratory wave rotor experiments. Both experiments 
utilize the so called divider cycle (Figure 4.1). In this 
cycle, the working gas (ah in both cases) enters the wave 
rotor through one port and exits through two others. In 
one of the exit ports the stagnation pressure of the gas is 
higher than that which enters the wave rotor, in the other 
it is lower. Although the practical utility of such cycles is 
debatable, they serve as excellent vehicles for wave rotor 
studies because they are relatively simple. In particular, 
they do not require external heating devices such as a 
burner, which would be needed in a full wave rotor gas 
generator experiment. Furthermore, wave rotor divider 
cycles clearly illustrate the concept of work exchange 
through unsteady wave propagation. 



Figure 4.1 Proposed Experimental Wave Diagrams 

The first experiment to be discussed is currently being 
conducted at NASA's Lewis Research Center. This wave 
rotor is highly instrumented making it possible to obtain 
both steady data in the ports and unsteady data from 
several x/L positions on one of the tubes as it rotates 
through the cycle. The data which will be presented was 
obtained from a single operating point. Unfortunately, just 
after the experiment began operation, the supply air to the 
entire facility was turned off for repairs. The experiment 
will resume in several months. The second experiment 
was performed by J. A. C. Kentfield in the mid 1960V 9> . 


There is no unsteady data available from this experiment; 
however, the steady results are extensive and as such, can 
be used to test the capabilities of the model over a broad 
range of conditions. Figure 4.1 shows the simplified wave 
diagrams which the two experiments were intended to 
follow. Table 4.1 shows the relevant geometrical 
dimensions and the gas conditions used in the 
experiments. It should be noted that the two experiments 
represent relative extremes of the various effects that are 
being modelled. Kentfield's passages had large cross 
sections, were short, and had small end clearances. This 
made the frictional and leakage effects small and the finite 
opening time effects large. The NASA experiment on the 
other hand, has passages with very small cross sections 
and relatively large leakage paths. 



Kentfield 

NASA 

Mean Tube Width, b 

0.60" 

0.25” 



0.40" 

Rotor O.D., D 

8.00° 

12.00" 

Rotor Length, L 

ii.oo" ; 

18.00" 

End Clearance, 5 

0.007" 

0.01" ports 
0.02" walls 

Rotor Speed, oo 


4150 rpm 


3 

i 

^om 

555 R 

637 R 

P 0M 



Pol 

14.70 psia 


a Equation 3.1.4 

2.579 

2.553 j 


Table 4.1 Experimental Parameters 

The model was tested in the following manner. The 
geometric parameters were entered into the code along 
with the inlet stagnation conditions and exit static 
pressures. Also entered was a guess at the pressure and 
temperature in the cavity of the wave rotor. These values 
are used in the leakage calculations. The code, which 
actually follows a single tube, was then run through a 
cycle. As it computed, a running sum of the mass 
momentum and energy fluxes in the ports at each time 
step was kept as well as a sum of the total leakage mass 
and energy flow into the cavity. At the end of the cycle, 
the cavity pressure and temperature were updated using 
the leakage sums in a simple lumped parameter model of 
the cavity region, and the cycle was run again. This 
process was repeated until the cavity pressure and 
temperature reached a steady value, and the total mass 
flow into the wave rotor matched the total outflow. At 
this point the integrated port fluxes were used to calculate 
a representative uniform fluid state, averaged using a 
constant area mixing calculation 0 *. This state was then 
used to obtain the outflow stagnation conditions. The 
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condition of steady, matched flow in the ports is termed 
the limit cycle of the wave rotor for a specified set of port 
boundary conditions. It is the wave cycle which any tube 
on the rotor will continually repeat. In comparing the 
model with the Kentfleld experiment a rather crude grid 
spacing of Ax/L=0.02 was used in order to reach the limit 
cycle condition as quickly as possible. The value of a in 
equation 3. 1.4 was chosen so as to obtain the best overall 
match to the data. For comparison with the present 
NASA experiment, a finer grid spacing of Ax/L=0.()l was 
used and the value of a was chosen such that the 
computed ratio of high pressure exit mass flow to input 
mass flow, p matched the experiment. The two values of 
a are listed for comparison in Table 4.1. it is seen that 
they are quite close, despite the very different geometries 
of the two experiments. This result lends some confidence 
to the analysis of the viscous effects. 



Figure 4.2 Unsteady Pressure Trace at x/L=.975 

Figure 4.2 shows a pressure trace measured in the tube of 
the NASA experiment, as it rotated through the cycle, at 
a location corresponding to 2.5% from the right side of 
Figure 4.1(b). All of the data has been non- 

dimensional ized by the inlet stagnation pressure. The 
circular symbols represent the experimental data, the solid 
line was obtained using the model with friction and 
leakage effects. The shaded regions indicate those 

portions of the cycle where the end of the tube nearest the 
transducer was adjacent to a wall. It can be seen that 
there is excellent agreement between predicted and actual 
results. Furthermore, it was found that the predicted mass 
flow through the rotor was only 8.3% higher than the 
measured How rate. This discrepancy might be accounted 
for by considering the blockage caused by the thickness of 
the tube walls. For reference, the dotted line in the figure 
represents the model output with only frictional effects. 
Although the timing of the waves appears correct the 
overall match with the experimental data is much worse. 
Furthermore, it was only possible to obtain the 
experimental mass flow split, p if the friction coefficient 
a was increased 2.24 times over the case with leakage 
modelled. This result indicates that leakage flows have a 
considerable effect on the wave pattern. 



Figure 4.3 Unsteady Pressure Trace at x/L=.250 

Figures 4.3 and 4.4 show the same comparisons between 
model and experiment for the x/L locations of .025 and 
.50 respectively, as measured from the left side of Figure 
4.1(b). Again, the agreement is quite good. It is noted 
that the predicted pressure in the tube at the end of the 
cycle is slightly lower than the experimental pressure in all 
three figures. Also, the predicted pressure in Figure 4.3, 
just before the opening of the inlet port is too high. These 
discrepancies may be due to improperly specified 
clearances (values used in the model were obtained from 
specs, not measured) or errors resulting from the 
assumptions of leakage model itself. 



Figure 4.4 Unsteady Pressure 7 race at x/L=.500 

Figure 4.5 shows the results of Kentfield’s experiment 
compared to the model predictions. Here, the ratio of the 
high pressure exit stagnation pressure to the inlet (medium 
pressure) stagnation pressure is plotted as a function of die 
ratio of the low pressure exit stagnation pressure to the 
inlet stagnation pressure for various values of the mass 
flow ratio (3 (high pressure port mass flow to total mass 
flow). The agreement between predicted and experimental 
results is again, quite good. The largest discrepancies 
appear in the lower right portion of the plot. In this 
region, the mass flow, and corresponding port velocities 
are quite low. The circumferential component of velocity 
brought about by the rotation of the tubes is actually larger 
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than the axial component generated by the waves and it is 
probable that this has a substantial effect on the 
performance. It is interesting to note that the lines of data 
all rise to the left and then abruptly stop. Kentfield 
pointed out that beyond this point the flow rate could not 
be increased. The same phenomenon was found in the 
model, although at slightly higher pressure. 



Figure 4.5 Kentfield’s Divider Performance Curves 
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5.0 CONCLUSIONS 

A numerical wave rotor model has been described which 
successfully predicted the behavior of two experiments 
with quite different geometries. The model is relatively 
simple in that it is one-dimensional, but sufficiently 
detailed to include the effects of viscosity, leakage and 
finite opening time. It requires relatively little computing 
time on a standard VAX mainframe (4.13 seconds CPU 
time for one cycle with 50 cells on a VAX 9410) and is 
therefore quite suitable for preliminary design purposes as 
well as analysis of existing wave rotors. 
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